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The discovery of the accelerated expansion of the Universe has had a vast resonance on a number 
of physical disciplines. In recent years several viable modified gravity models have been proposed, 
which naturally lead to a late-time de Sitter stage while basically reducing to General Relativity 
in the early Universe. We consider a contracting cloud of pressureless dust, having arbitrary mass 
and initial density, and discuss some aspects of these modified gravity models. We show how 
the increasing energy /mass density may lead to a curvature singularity and discuss the typical 
timescales for its development. 
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I. INTRODUCTION 

The physical mechanism behind the present accelera- 
tion of the cosmological expansion [l[ is still unknown. 
A pure cosmological constant term is quite natural both 
in General Relativity and in Quantum Field Theory, but 
then there remains the challenge of explaining the present 
value of A, in particular its smallness and the fact that 
il m ~ f^A (coincidence problem). 

More complicated but popular scenarios involve 
"quintessence" [|| models, in which a scalar field coupled 
to gravity is responsible for the acceleration, or modified 
gravity models, in which the acceleration is due to mod- 
ifications of the Einstein-Hilbert action and therefore of 
the Einstein field equations. In the simplest case, the 
gravitational lagrangian is a non-linear function of the 
scalar curvature R alone0 
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In this paper, we discuss the contraction of a nearly- 
homogeneous cloud of pressureless dust in f(R) gravity 
models. The arguments presented here should be consid- 
ered an extension to those in 0] , where the appearance of 
singularities was shown numerically and discussed semi- 
analytically, and in [4|. Also, there is a strong relation 
between the singularities dealt with in this paper and 
those considered in past cosmological scenarios, e.g. in 
refs. [U, [|| . From a slightly different perspective, grav- 
itational contraction in f(R) gravity has been recently 
discussed in Q. 
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We use natural units 
defined as m? pl = 
and we use the conventions 

R < for a matter-dominated Universe. 



- h = k = 1, and the Planck mass is 
The metric has signature (+ — ), 

(fl/tf/h, + ■■■). «v = 

With these conventions, 
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We find that the increase in the energy /mass density 
may indeed lead to curvature singularities, and that these 
are unavoidable if the contraction lasts long enough. 
Most results are quite general, and can be applied to 
a few popular models |8l— FlQj| constructed to generate a 
late-time de Sitter stage. 

Usually, singularities are assumed to be prevented by 
the addition of ultraviolet corrections to the gravitational 
lagrangian [3|, |6j , for instance quadratic terms which are 
natural when one-loop quantum corrections are consid- 
ered . A fully non-perturbative approach, taking into 
account all possible higher-order corrections, may shed 
more light on singularities in gravitational theories, or 
lack thereof [HJ . 

Moreover, an oscillating scalar curvature gives rise to 
gravitational particle production, whose back-reaction on 
gravity is a damping of the oscillations of R [l4| . This can 
also help avoiding the singularity, and could in principle 
yield detectable cosmic ray signatures (l5| . 



II. CURVATURE EVOLUTION IN 
CONTRACTING SYSTEMS 

From eq. (HJ, one obtains the field equations 



f(R)R^--f(R)g^ + 



(g^V 2 -T>rV v ) f{R)=T tiV . 



(2) 



Here, T> denotes covariant derivative, a prime denotes 
derivative with respect to R, and 



where A m is the matter action. The corresponding trace 
equation reads 




W 2 F' + RF' - 2F - (R+ T) = . 



(3) 



We will consider a nearly homogeneous, spherically sym- 
metric cloud of pressureless dust, hence T — 8irg/m Pl . 
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It can also be shown that spatial derivatives can be ne- 
glected [3]. Moreover, we assume to be in low-gravity 
regime, so that \F/R\ <C 1, \F'\ <C 1, the spacetime 
is nearly Minkowski, and covariant derivatives reduce in 
first approximation to ordinary derivatives. Indeed, these 
arguments can also be applied to cosmology, if we con- 
sider the evolution of the Universe, backwards in time, 
during the matter-dominated epoch. 
Under these assumptions, eq. (J3)) becomes 



3<9 2 F' - (R + T) = . 
Defining the new scalar field 

i = -3F' , 
we rewrite it as an oscillator equation: 



i + (R + T) = o £ 



dU 



(4) 



(5) 



(6) 



In general, it is not possible to invert ((SJ) to obtain 
R = R(£) and thus a simple form for U(£), expect per- 
haps in some limit, but it is rather clear that solutions 
will oscillate around the GR solution R + T = 0, with 
frequency roughly given by 



d 2 U 



i 



R+T=0 



d^/dR 



(7) 



R+T=0 



If uj 1 < 0, one expects instabilities, and this is exactly 
the kind of instability first discovered by Dolgov and 
Kawasaki [l3|]. Even with w 2 > 0, we will show that 
if the model fulfils some requirements, then curvature 
singularities can be developed. In particular, we need 
that 

• there exist a certain value losing corresponding to 
\R\ -> oo, 

• the potential be finite in £ s , ns , that is U{^ S i ng ) < 
oo. 

If the previous requirements are met, then in general it 
is possible that £ reach £ S i„ g and hence \R\ — > oo. We 
can see this, for instance, from the "energy" conservation 
equation associated to flBJ, that is 



1 



dT 

e + U(0- I dt'^m= const. (8) 



The last term is due to the (possible) explicit time- 
dependence of T, and if dT/dt > 0, as is the case in con- 
tracting systems, it will in general produce an increase in 
the "canonical" energy (kinetic + potential). Note that 
this is true for £ > 0, whereas for £ < it would give the 
opposite behaviour. However, it has been shown that the 
condition F' < 0, corresponding to £ > 0, is crucial for 
the correct behaviour of the modified gravity model at 
(relatively) low curvatures [6j. 



For simplicity, let us assume that the oscillations of £ 
in its potential are "adiabatic" , in the sense that at each 
oscillation £ moves between two values ^ m in and £ max 
at roughly the same "height", that is with U(£ m i n ) — 
U(£ max ). Then, if at some moment U(£) = U(^ sing ), 
we expect that £ — ► £ S i„ ff within at most the follow- 
ing oscillation. In what follows, we will show that this 
is a rather natural result, provided that the contraction 
lasts long enough, or better that the contraction factor 
Tfinai/Tinitiai is large enough. 

We expand £ around the "average" value £ = £ Q (i) 
corresponding to the GR solution, namely 



where 



and 



= &(t) + a(t) sin $(*), 



i?(£ Q )+T = 



$(t) ~ / dt'u 



(9a) 
(9b) 



(10) 



(11) 



The function a is assumed to be relatively slowly- varying, 
that is a/a < uj. 



A. First-Order Expansion 

At first order in £i (defined in eg. I9"aj) . eq. (JBJ) reads 

& + w 2 & (12) 

Using the expansion (|9b[) and neglecting £ a and a yields 



uj a 



a(t) ~ ao 



w(t) 



(13) 



As long as the approximations hold, this can be con- 
sidered a rather general result, and the specific F(R) 
model will determine the behaviour of the oscillations. 
The value ao in eq. (fT3")) is strictly related to the initial 
conditions we impose, that is to the initial displacement 
from the GR behaviour. We will in general set the initial 
values of R and R, and from those derive the initial val- 
ues of £ and £. Thus, ao can be calculated differentiating 
eq. I|9bp. which yields 

£o(Rq, Rq) ^ £a,o + «o^o => «o - — — — ■ (14) 

This corresponds to the explicit solution 

"(i)-((o-ea,o)hw(tf 1/2 . (15) 

Please note that, apparently, we have not made use of the 
assumption U(^ s i ng ) < oo. Although not necessary to 
perform calculations, this condition is needed to ensure 
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that expansion ([9]) be reliable. In fact, oscillations are 
harmonic only if the potential is nearly quadratic; this 
assumption is usually quite reasonable, especially near 
the bottom of the potential, but loses its validity, for 
instance, near points at which U diverges. Therefore, 
models in which 



(16) 



is singular in £ = ( si „ 9 , unlike the three models [§4l0| we 
will explicitly consider later, cannot be discussed within 
the framework of this paper. 

Also, it is clear from eq. (IT51) that if £ = £ Q the am- 
plitude of oscillations should vanish at all times. This 
can be immediately proved wrong, for instance numeri- 
cally. The disagreement with our estimate is due to the 
fact that we had neglected £ Q , and in general terms pro- 
portional to l/t 2 contrl in eq. (H"2j) . When a is initially very 
small, however, those terms should be kept and the ap- 
proximations used are no longer valid. Therefore, eq. (|15[) 
is reliable when (£ — £ a ,o) is "large" enough, say of the 
order of £ a ,o- 



III. GENERATION OF CURVATURE 
SINGULARITIES 

The singularity develops as soon as £ = £, S ing, or soon 
after [/(£) = U(£ s i ng ). Exploiting the results of the pre- 
vious paragraph, we see that this is basically equivalent 
to the condition 



\Ut)-Ung\ <a{t) 



(17) 



Therefore, the singularity should be reached roughly at 
the moment when condition (I17[) is fulfilled. Substituting 
(fl3l) results in 



\ut)-i 



sing 



< 



(wwo) 



-1/2 



(18) 



In order to discuss these results more quantitatively, we 
shall apply them to three models which have recently 
been proposed: 



F(R) 



XRr 



l + (i?/i? c )- 2 "' 
)2\ -™ 



F(R) = XR C 
F(R) = R c ]n 
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IF 
Rl 



(19a) 
(19b) 



- A + (1 - e- x ) e- R / R < 



m ( 19c ) 



Note that for all three models, if A is of order unity R c is 
of the order of the present cosmological constant, which 
is much smaller than the typical values of R and T in 
astrophysical systems, such as pre-stellar, pre-galactic, 
and molecular clouds. Hence, in many cases we will take 
the limit \R C /T\ ~ R c /R <C 1 before presenting the final 
results. 



For simplicity, we assume that the contraction of the 
system is stationary, i.e. the mass density grows linearly 
with time: 



T(t) = T (1 + t/t contr ) 



(20) 



This evolution law is not completely accurate from a 
physical standpoint, but unless the contraction follows a 
very different behaviour, and in any case until t < t con t r , 
which is assumed to be the typical timescale of contrac- 
tion, results obtained with this form should be more or 
less accurate. In order to have simple and somewhat re- 
liable estimates, we will use the initial conditions 



(21) 



Depending on physical and model parameters, we will 
calculate the "critical" energy density T S i ng , at which 
the singularity arises. The timescale for the formation of 
the singularity, as can be seen from eq. (1201) , is in general 




t s 



T s „ 



To 



(22) 



As previously stressed, the evolution law (|20|) can be con- 
sidered somewhat reliable until t is of the order of magni- 
tude of t contr , so we will focus on combinations of param- 
eters in which T S i ng /To is at most a few. Furthermore, 
we will define the following quantities: 



R. 



-29 = -J — 



\R C 



8tt 10" 29 g cm" 3 



£•29 



10~ 29 g cm~ 3 



(23) 
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J !0 = TTtIO - 
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A. Hu-Sawicki and Starobinsky 

We will show explicit expressions for the Hu-Sawicki 
model [K|. This model and that by Starobinsky || are 
essentially identical in the limit \R\ ^> \R C \, so the results 
of this paragraph can be applied to both models. The 
field £ is 



£(R) = 6nA 



(R/Rcf 



[i + (R/R c y 



(24) 



therefore \R\ — > oo corresponds to £ s i ng — 0. We also 
have 



Rc 



T 



2n+2 



6n\(2n + 1) \ R c 
so with the initial conditions (|2"Tj) we have 

j, x -(2n+l) 



Ca - 6nA ( - — 



(25) 



(26) 
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Ua _**^rg\-™ (27) 

tcontr V / 



'1/2 



/6nA(2n+ 1) /TT 



\Rc. 



-(n+l)/2 



(28) 



Substituting these expressions and using the definitions 
(pi), eq. (p]l yields 



X Q L contr 



2/(3n+l) 



o V V 6 « A ( 2 " + 1 ) 3i? " +1/2 , 



0.53- 



029 HO 



2/(3n+l) 



V"A(271 + l) 3 A 



n+1/2 
29 



(29) 



Typically we would have A, R29 ~ C(l), so from eq. ([29]) 
we see that if the initial energy density is not too large 
and/or the contraction is fast enough, the overall con- 
traction needed to reach the singularity is quite realis- 
tic, with T sing /T possibly even only slightly larger than 
unity. Moreover, the larger n, the easier it is to develop 
the singularity with relatively small T S i ng . 



B. Appleby-Battye 

In this case the field £ is 

We also find 



(30) 



R c 



3(e A - 1) 



-T/R c 



£.a,0 



3(e A -I) ( T ( 



tcontr 

Therefore, ([18$ yields 
5 T AR, 



(31) 



T, 



sing 



In - 



■f) 11 



V3(e A - 1)T 

r3/2 



(32) 



which leads to 



T s! 



To 



f_R2 

3 £>29 



In 



0.9We x - lg 



2<> 



p3/2 , 



(33) 



As in the previous models, reasonable values of physical 
and model parameters may very well lead to a singular- 
ity after T has grown by a small factor from its initial 
value. Indeed, due to the weak logarithmic dependence 



and since usually one assumes To ^> 
very close to 5Tq/3. 



Rc\, T sing is usually 



IV. CONCLUSIONS 

The possibility of curvature singularities in F(R) grav- 
ity models in the presence of contracting bodies has been 
discussed in a rather simple fashion. The trace of the 
modified Einstein equations has been recast as an oscil- 
lator equation for the field £ ~ dF/dR, which moves in a 
time-dependent potential U (£, t). This equation has been 
solved expanding it at first order in deviations from GR. 
Generally speaking, solutions oscillate around the "aver- 
age" value, corresponding to the GR solution R + T = 0, 
with increasing amplitude. 

Models in which \R\ — > 00 implies neither £ — > 00 nor 
U — > 00 are particularly suited to be discussed within this 
framework, and this is the case for three well-known mod- 
els (8l-[To|. In all three cases, a curvature singularity is 
to be expected if the mass/energy density reaches T s i ng , 
which depends both on model parameters and on physical 
properties of the system under consideration. For typical 
values of model parameters and reasonable physical sys- 
tems, the singularity is expected on timescales possibly 
much shorter than the age of the Universe. 

Two mechanisms should help avoiding the singular- 
ity, namely ultraviolet gravity modifications and grav- 
itational particle production. The results of this work 
could be used both to estimate the cosmic ray flux due 
to gravitational particle production, and possibly to con- 
strain modified gravity models at large curvatures [l5| . 
This could be the subject of future research. 
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